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ABSTRACT 

The survey of galaxy clusters performed by Planck through the Sunyaev-Zeldovich effect has already discovered many interesting objects, thanks 
to the whole coverage of the sky. One of the SZ candidates detected in the early months of the mission near to the signal to noise threshold, 
PLCKG214.6+37.0, was later revealed by XMM-Newton to be a triple system of galaxy clusters. We have further investigated this puzzling system 
with a multi-wavelength approach and we present here the results from a deep XMM-Newton re-observation. The characterisation of the physical 
properties of the three components has allowed us to build a template model to extract the total SZ signal of this system with Planck data. We 
partly reconciled the discrepancy between the expected SZ signal from X-rays and the observed one, which are now consistent at less than 1.2 cr. 
We measured the redshift of the three components with the iron lines in the X-ray spectrum, and confirmed that the three clumps are likely part 
of the same supercluster structure. The analysis of the dynamical state of the three components, as well as the absence of detectable excess X-ray 
emission, suggest that we are witnessing the formation of a massive cluster at an early phase of interaction. 

Key words. Cosmology - Clusters of galaxies 



1. Introduction 

Clusters of galaxies occupy a special position in the hierarchy 
of cosmic structures: they are the largest objects that decoupled 
from the cosmic expansion and that have had time to undergo 
gravitational collapse. They are thought to form via a hierarchi- 
cal sequence of mergers and accretion of smaller systems driven 
by gravity. During this process the intergalactic gas is heated to 
high X-ray emitting temperatures by adiabatic compression and 
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shocks and settles in hydrostatic equilibrium within the cluster 
potential well. Sometimes galaxy clusters are found in multiple 
systems, super-cluster structures which already decoupled from 
the Hubble flow and are destined to collapse. The crowded envi- 
ronment of superclusters is an ideal place to study the merging 
processes of individual components at an early stage of merger 
and witness the initial formation phase of very massive struc- 
tures. Moreover, the processes related to the contraction may in- 
crease the density of the intercluster medium and make it ob- 
servable with present instruments. An example is the central 
complex of the Shapley concentration which has been the ob- 
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ject of extensive multi-wavelength observations with the aim 
of characterising the merger processes in galaxy clusters (Kull 
^^5Mnger[{T999 l Bar delli et alT1[T998) |Rossetti et al.||2007| 
Giacintucci et al. 20Q5|>. 

Recently a new observational window has opened up for the 
study of the astrophysics of galaxy clusters through the Sunyaev- 
Zeldovich effect ( |Sunyaev & Zeldovich |1972| SZ hereafter): 
a spectral distortion of the Cosmic Microwave Background 
(CMB) generated through inverse Compton scattering of CMB 
photons by thermal electrons in the intracluster medium (ICM). 
SZ surveys are discovering new clusters, some of which are in- 
teresting merging systems (with "El Gordo", Menanteau et al. 
|201 1[ being probably the most spectacular example). 
A key role in SZ science is now played by the PZohcAQ satel- 
lite. Compared to other SZ surveys of galaxy clusters, Planck 
has only moderate band-dependent spatial resolution, but it pos- 
sesses a unique nine band coverage and more crucially it covers 
the whole sky. Therefore, it allows the detection of the rarest 



2. Observations 



2. 1 . Planck data and analysis 



objects: massive high-redshift systems (Planck Collaboration 
201 If) , which are the most sensitive to cosmology, and com- 



plex multiple systems, which are interesting for the physics of 
structure formation. Indeed, during the follow-up XMM-Newton 
campaign of Planck SZ candidates, we found two new dou- 



ble systems and two new triple systems of clusters (Planck 



Collabora tion|201 lc hereafter Paper I). In all cases, the cumu- 
lative contribution predicted by X-ray measurements was lower 
than the measured SZ signal, although compatible within three 
cr. 

PLCKG2 14.6+37.0 is the most massive and the X-ray brightest 
of the two Planck discovered triple systems. The XMM-Newton 
follow-up observations showed that the Planck SZ source candi- 
date position is located ~ 5' from the two southern components 
(A and E). A third subcomponent, C, lies approximately 7' to 
the North (Fig. [TJ. The X-ray spectral analysis of the compo- 
nent A indicated a redshift of zpe ~ 0.45, consistent with two 
galaxies with spectroscopic redshift of ~ 0.45, close to the peaks 
of component A and C. A cross-correlation with SDSS-DR7 
Luminous Red Galaxies and the Superclusters catalogue from 
the SDSS-DR7 ( |Liivamagi et al.|2010| > hinted that this triple sys- 
tem is encompassed within a very large-scale structure located 
at z ~ 0.45, and whose centroid lies about 2° to the South (see 
Appendix B in Paper I for further details). 
In this paper we present new SZ measurements of this object 
with Planck and compare them with the results from a a deep 
XMM-Newton re-observation. In Sect. 2 we describe the analy- 
sis methods and the Planck and XMM-Newton data used in this 
paper. In Sect. 3 we present our results obtained with X-ray and, 
in Sect. 4 compare them with available optical data from SDSS. 
In Sect. 5 we compare the X-ray results with Planck results. In 
Sect. 6 we discuss our findings. 

Throughout the paper we adopt a ACDM cosmology with Ho = 
70 km s" 1 Mpc" 1 , Q M = 0.3 and Q A = 0.7. At the nominal 
redshift of the supercluster, z = 0.45, one arcminute corresponds 
to 347 kpc. 



1 Planck (http://www.esa.int/Planckl is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries: France and Italy) with contributions from NASA (USA), and 
telescope reflectors provided in a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



Planck (Taube rlst~al.||2010[ |Planck Col laboratio n|201 la| > is the 
third-generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI ; JMandolesi et al.pOTOt |Bersanelli et~aT1[20T0l 
|Mennella et al.||2011| i covers the 30, 44, and 70 GHz bands 
with amplifiers cooled to 20 K. The High Frequency Instrument 
(HFI; |Lamarre et al.|2010]|Planck HFI Core Team|201 la] > covers 
the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0. 1 K. Polarisation is measured in all but the 
highest two bands ( |Leahy et aLl|2010| |Rosset et aLl|2010| >. A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics ( |Planck Co llaborati on|201 lb) . Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
( |Planck HFI Core Team||2011b| |Zacchei et a"Ll[20TT] i. Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration, 201 lh-z. 

Our results are based on the SZ signal as extracted from the six 
bands of HFI corresponding to the nominal Planck survey: 14 
months, during which the whole sky was observed twice. We re- 
fer to |Planck HFI Core Team] ( [20 lib) and |Zacchei eTaE] ( |20TT) 
for the generic scheme of TOI processing and map making, as 
well as for the technical characteristics of the maps used. We 
adopted a circular Gaussian as the beam pattern for each fre- 
quency as described in |Planck HFI Core Team (201 lb}; |Zacchei| 
|etlT1 ( |20TT) . 

The total SZ signal is characterised by the integrated Compton 
parameter Y500 defined as D 2 A (z)Ysao = (cr T lm e c 2 ) J P,i,dV, 
where Da(z) is the angular distance to a system at redshift z, 
o~t is the Thomson cross-section, c the speed of light, m e the 
rest mass of the electron, P,/, is the pressure of thermal electrons 
and the integral is performed over a sphere of radius P500. 
The extraction of the total SZ signal for this structure is more 
complicated than for single clusters. Due to its moderate spa- 
tial resolution, Planck is not able to separate the contributions 
of the three components from the whole signal. In Paper I, we 
estimated the total flux assuming a single component with mass 
corresponding to the sum of the masses of the three clumps, fol- 
lowing a universal pressure profile ( Arnaud et al.|2010| centred 



at the barycentre of the three components. This is obviously a 
simple first-order approach. The wealth of information that is 
now available on this system, has allowed us to build a more 
representative model. With the X-ray constraints on the struc- 
tural properties of clumps A, B and C (Sec. |3.1[ ), we are now 
able to build a three-component model. As discussed in Planck 
[Collaboration ( 201 ld| l, our baseline pressure profile is the stan- 
dard "universal" pressure profile derived by Arn aud et al.| ( |2010] l. 
We assumed this profile in each clump, parametrized in size by 
the respective X-ray scale radius, Rsm- The normalisations, ex- 
pressed as integrated Comptonization parameters within 5/?5oo 
were tied up together according to the ratio of their respective 
ix,500, as determined from the X-ray analysis. Thus only one 
overall normalisation parameter remains to be determined. This 
template used under these assumptions and parametrisations to- 
gether with the multi-frequency matched filter algorithm, MMF3 
( |Melin et al.|200 6 ), directly provides the integrated SZ flux over 
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Fig. 1: The triple systems PLCKG2 14.6+37.0. Left panel: Planck SZ reconstructed map (derived from the Modified Internal Linear 
Combination Analysis, Hurier et al. |2010 l, oversampled and smoothed for displaying purposes. Middle panel: XMM-Newton 
wavelet filtered image in the [0.5-2.5] keV (Sec 2.2). The three components of PLCKG214. 6+37.0 are marked with circles (with 
radius R500, Table 1) and indicated with letters as in the text. The X-ray structure visible on the left (marked with *) is likely associ- 
ated with the active galaxy SDSSJ090912.15 + 144613.7, with a spectroscopic redshift z = 0.767. Right panel: SDSS r-band image 
of the PLCKG2 14.6+37.0 field. In all panels, North is up and East to the left and X-ray contours are overlaid. 



the whole super-structure. 

The two-dimensional reconstruction of the Comptonization pa- 
rameter, y = (crj/meC 2 ) f Pfkdl, provides a way to map the spa- 
tial distribution of the thermal pressure integrated along the line 
of sight. This is performed using the modified internal linear 
combination method (MILCA, |Hurier et al.||2010] > on the six 
Planck all-sky maps from 100 GHz to 857 GHz. MILCA is a 
component separation approach aimed at extracting a chosen 
component (here the thermal SZ signal) from a multi-channel 
set of input maps. It is mainly based on the ILC approach (e.g., 
|Eriksen et al.|2004) , which searches for the linear combination 
of the input maps that minimises the variance of the final recon- 
structed map imposing spectral constraints. 

2.2. XMM-Newton observation and data reduction 

PLCKG2 14.6+37.0 was re-observed by XMM-Newton during 
AO10, for a nominal exposure time of 65 ks. We produced cal- 
ibrated event files from Observation Data Files (ODF) using 
v.11.0 of the XMM-Newton Science Analysis Software (SAS). 
We cleaned the event files from soft protons flares, using a dou- 
ble filtering process (see Bourdin & Mazzotta|200 8 for details). 
After the cleaning, the net exposure time is ~ 47 ks for the MOS 
detectors and ~ 37 ks for the pn. Since a quiescent component 
of soft protons may survive the procedure described above, we 
calculated the "in over out ratio" 7?sb ( |De Luca"& Molendi 2004 ) 
and we found values close to unity, suggesting a negligible con- 
tribution of this component. We masked bright point sources, de- 
tected in the MOS images as described in Bourdin & Mazzotta 
(2008 ). We performed the same data reduction procedure on the 
snapshot observation 0656200101, finding a cleaned exposure 
time of ~ 15 ks for the MOS and ~ 10 for the pn. 
We then combined the two observations and binned the pho- 
ton events in sky coordinates and energy cubes, matching the 
angular and spectral resolution of each focal instrument. For 
spectroscopic and imaging purposes, we associated an "effec- 



tive exposure" and a "background noise" cube to this photon 
cube (see Bou rdin et al.|[20TT| for details). The "effective ex- 
posure" is computed as a linear combination of CCD exposure 
times related to individual observations, with local corrections 
for useful CCD areas, RGS transmissions, and mirror vignetting 
factors. The "background noise" includes a set of particle back- 
ground spectra modelled from observations performed with the 
closed filter. Follo wi ng an approach proposed in e. g. , |Leccardi| 
& Molendi| ( |2"0"0"8"l > or |Kuntz & Snowden| ( (2008l l, this model sums 
a quiescent continuum to a set of fluorescence emission lines 
convolved with the energy response of each detector. Secondary 
background noise components include the cosmic X-ray back- 
ground and galactic foregrounds. The cosmic X-ray background 
is modelled with an absorbed power law of index F = 1 .42 (e. g. , 
Lumb et al.|2002 1, while the galactic foregrounds are modelled 



by the sum of two absorbed thermal components accounting 
for the galactic transabsorption emission ( kT\ — 0.099 keV, 
kT 2 = 0.248 keV, |Kuntz & Snowden|2000| l. We estimate emis- 
sivities of each of these components from a joint fit of all back- 
ground noise components in a region of the field of view located 
beyond the supercluster boundary. 

To estimate average ICM temperatures, kT along the line of sight 
and for a given region of the field of view, we added a source 
emission spectrum to the "background noise", and fitted the 
spectral shape of the resulting function to the photon energy dis- 
tribution registered in the 0.3 - 10 keV energy band. In this mod- 
elling, the source emission spectrum assumes a redshifted and 
nH absorbed emission modelled from the Astrophysical Plasma 
Emission Cod e (APEC, |Smith et al.||2001| l, with the element 
abundances of Grevesse & Sauval|(|1998|l and neutral hydrogen 



absorption cross sections of |Balucinska-Church & McCammon 



( |1992| l. It is corrected for effective exposure, altered by the mir- 
ror effective areas, filter transmissions and detector quantum ef- 
ficiency, and convolved by a local energy response matrix. 
The X-ray image shown in Fig. |T| is a wavelet filtered image, 
computed in the 0.5 - 2.5 keV energy band. To generate this 
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image, we corrected the photon map for effective exposure and 
soft thresholded its undecimated B3-spline wavelet coefficients 
( Starck et al. 2001) to a 3 cr level. In this procedure, significance 
thresholds have been directly computed from the raw (Poisson 
distributed) photon map, following the multi-scale variance sta- 
bilisation scheme introduced in Zhang et al. ( 2008| l. We applied 



the same transformation to a "background noise" map, which we 
then subtracted from the image. 



3. Structure of the clusters from X-rays 

3. 1 . Global analysis of the cluster components 

Assuming that the three structures are located at the same red- 
shift, z = 0.45 (the spectroscopic value found in Paper I), from 
the combination of the two XMM-Newton observations we have 
carried out an X-ray analysis for each component independently 
(we masked the two other clumps while analysing the third 
one). We extracted surface brightness profiles of each compo- 
nent, centred on the X-ray peak, in the energy band 0.5 - 2.5 
keV. We used the surface brightness profile to model the three- 
dimensional density profile: the parametric density distribution 
( |Vikhlinin et al.|[2006[ ) was projected, convolved with the PSF 
and fitted to the observed surface brightness profile (Bourdin 
et al. 201 1] >. From the density profile we measured the gas mass, 
M g , that we combined with the global temperature Tx, obtained 
with the spectral analysis, to measure Y x = M g * T x ( |Kravtsov 
et al.||2006|l. We used the Msoo-Yx scaling relation in |ArnaucT 
et al-| ( |2010| > to estimate the total mass M500, defined as the 
mass corresponding to a density contrast 6 = 500 with respect 
to the critical density at the redshift of the cluster, p c (z), thus 



M500 = (47r/3)500p c (z)/?j 00 . The global cluster parameters were 
estimated iteratively within R500, until convergence. 
The resulting global X-ray properties are summarised in Table[T] 
The Yx and M500 values are slightly larger, but consistent within 
less than 1.5 cr with the results shown in Paper I. 



3.2. Redshift estimates 

Crucial information on the nature of this triple system comes 
from the measurement of the redshift of each component, allow- 
ing us to assess whether this is a bound supercluster structure or a 
combination of unrelated objects along the same line of sight. In 
Paper I, a reliable redshift measurement, obtained with the short 
XMM-Newton observation, was available only for component 
A. Its value (z = 0.45) was consistent with the only two spec- 
troscopic redshifts available in this field and corresponding to 
the bright central galaxies in components A (SDSSJ090849.38 + 
143830.1 z = 0.450) and C (SDSSJ090851.2 + 144551.0 z = 
0.452). A photometric redshift, z = 0.46, was furthermore avail- 
able for a bright galaxy (SDSS J090902.66+ 143948. 1) very close 
to the peak of component B. 

With the new XMM-Newton observation, we detect the iron K 
complex in each clump. We extracted spectra in a circle cen- 
tered on each component with radius R50Q (Table [TJ. We per- 
formed a more standard spectral extraction and analysis than the 
one described in Sect. 



2.2 extracting in each region the spec- 



trum for each detector and its appropriate response (RMF) and 
ancillary (ARF) files. We fitted spectra within XSPEC, mod- 



elling the instrumental and cosmic background as in Leccardi & 
Molendi (2008), leaving as free parameters of the fit the temper- 
ature, metal abundance, redshift and normalisation of the cluster 
component (see Planck Collaboration 201 lc| for details). We first 
fitted spectra for each detector separately. While the MOS detec- 



tors do not show any instrumental line in the whole 4-5 keV 
range (Leccardi & Molendi 2008), the pn detector shows a faint 
fluorescent lin^] in the spectral range where we expect to find 
the redshifted cluster line. We verified that this feature does not 
affect significantly our results, since the pn redshift and metal 
abundance are always consistent with at least one of the MOS 
detectors. We report our results in Table [2] 
For components A and C the redshift measurements for each de- 
tector are consistent within one <x and we performed joint fits 
combining all instruments (Table pi. For component B, the red- 
shift measurement with MOS 1 is larger than the estimates with 
the other detectors, although consistent within two cr. The joint 
fit of the three detectors in this case would lead to a best fit value 
z = 0.5 16 (-0.023, +0.014), while combining only MOS2 and 
pn we find 0.481 (-0.011, +0.013). In the joint fit of the three 
detectors, the MOS 1 spectrum drives the redshift estimate, leav- 
ing many residuals around the position of the iron lines for the 
MOS2 and pn spectra. We performed simulations within XSPEC 
to quantify the probability that, given the statistical quality of 
the spectra and the source to background ratio, a redshift mea- 
surement as large as z = 0.529, may result just from statistical 
fluctuations of a spectrum with z = 0.48. We assumed the best 
fit model of the joint MOS2 and pn analysis as the input source 
and background model with a redshift z = 0.48 and we gener- 
ated 1500 mock spectra that we fitted separately with the same 
procedure we used for real spectra. We found a redshift as large 
as what we measured with the MOS1 in 3% of the simulated 
spectra. With these simulations, we reproduced also the joint fit 
procedure: we performed 500 joint fits of three simulated spec- 
tra and found that a redshift as large as 0.516 occurs with less 
than 1% probability. Furthermore, we performed other simula- 
tions assuming the redshift resulting from the joint fit of the three 
detectors z = 0.516: the probability of finding two redshifts as 
low as z = 0.475 is about 1.7%. The simulations we have de- 
scribed show that it is unlikely that the three redshift measure- 
ments and the joint fit we obtained for subcluster B may result 
from statistical fluctuations of the same input spectrum, suggest- 
ing a systematic origin for the discrepancy between MOS 1 and 
the other detectors. Indeed, the quality of the fit with the MOS1 
data alone is worse than with the other detectors, featuring many 
residuals around the best fit model. We verified the possibility 
of a calibration issue affecting MOS1 by checking the position 
of the instrumental lines: we did not find any significant sys- 
tematic offset for the bright low-energy Al and Si lines but the 
absence of strong fluorescent lines between 2 and 5.4 keV did 
not allow us to test the calibration in the energy range we are 
interested in. Although the origin of the systematic difference of 
the MOS 1 spectrum is still unclear, we decided to exclude this 
detector when estimating the redshift of component B (Table [2 1. 
Nonetheless, in the following, we will also discuss the possibi - 
ity that the cluster is located at the larger redshift z = 0.516. 
Concerning the components A and C the redshift measurements 
are not significantly affected if we exclude the MOS 1 detector. 
The redshift estimates we obtained from X-ray data for compo- 
nents A (z = 0.445 + 0.006) and C (z = 0.46 + 0.01) are nicely 
consistent with the spectroscopic values found in the SDSS 
archive for their central brightest galaxies (0.450 and 0.452, 
respectively). Concerning component B, even without consid- 
ering the MOS1 detector, we still find a larger best fit value 
(z = 0.48 ± 0.01) with respect to the other two components. 
This is shown in Fig. [2J where we compare the variation of x 2 



2 Ti K„, http://xmm2.esac.esa.int/external/ximti_sw_cal/ 
background/f ilter_closed/pn/m£reyberg-WA2- 7 . ps . gz 
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Table 1: Physical properties of the components of PLCKG2 14.6+37.0 . 



Component 


RA X 


Dec x 




M e ,500 




M,oo 


-^500 




[hh:mm:ss] 


[hh:mm:ss] 


[keV] 


[10 14 M o ] 


[10 14 M o keV] 


[10 14 M o ] 


[kpc] 


A 


09 : 08 : 49.6 


+ 14 : 38 : 26.8 


3.6 + 0.4 


0.26 ± 0.01 


0.96 + 0.11 


2.22 + 0.16 


784+ 19 


B 


09 : 09 : 01.8 


+ 14 : 39 : 45.6 


4.3 + 0.9 


0.28 ± 0.02 


1.2 + 0.3 


2.5 ± 0.4 


820 ± 40 


C 


09 : 08 : 51.2 


+ 14 : 45 : 46.7 


5.3 + 0.9 


0.30 ± 0.02 


1.6 + 0.3 


3.0 + 0.3 


864 ± 33 



Table 2: Redshift measurements from the X-ray iron line for the components of PLCKG2 14.6+37.0. 



Component MOS1 MOS2 pn joint fit 



A 0.447 (-0.013 +0.024) 0.446 (-0.005 +0.013) 0.441 (-0.009 +0.009) 0.445 (-0.006 + 0.006) (ml +m2+pn) 
B 0.529 (-0.018 +0.024) 0.472 (-0.008 +0.063) 0.475 (-0.016 +0.023) 0.481 (-0.011 + 0.013) (m2+pn) 
C 0.469 (-0.020 +0.020) 0.434 (-0.016 +0.017) 0.463 (-0.017 +0.009) 0.459 (-0.010 + 0.010)(ml+m2+pn) 
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Fig. 2: Variation of x 1 when fitting the spectra for the redshift 
measurement. Black, red, and green lines represent component 
A, B, and C, respectively. Dashed thin lines correspond to the 
68% and 90 % error range, respectively. 



for the joint fits (see Table |2]i for the three components. While 
component A and C are consistent with being at the same red- 
shift at less than one <x, component B is likely located at a larger 
redshift, although consistent at less than two cr with the posi- 
tion of the other clusters. Therefore, component B is likely sep- 
arated along the line of sight from the two other components by 
69 (-30, +25) Mpc (150 Mpc, if we consider the redshift esti- 
mate obtained with the three detectors). While this large separa- 
tion suggests that the cluster B is not interacting with the other 
components, it is still consistent with the three objects being part 
of the same supercluster structure (Bah call|1999[ l. 
Using the best fit redshift estimates for the three components, we 
recomputed the physical parameters in Table 1 . While the vari- 
ations for cluster A and C are negligible, for cluster B we found 
Y x = (0.74 ± 0.19) 10 14 Mo keV, M 500 = (1.89 ± 0.20) 10 14 M o 
and Rsqo = (726 + 25) kpc. 

3.3. Radial structural analysis of the components 

We performed a radial analysis of the X-ray observations of 
PLCKG2 14.6+37.0 to study the behaviour of the main thermo- 
dynamical quantities of the ICM. We extracted the surface 



brightness profile as discussed in Sec. |3.1| while the A com- 
ponent shows a very peaked profile, which might indicate a 
cool core state, the B and C components have flatter profiles 
at the centre, a signature of an un-relaxed dynamical state. On 
more quantitative bases, we extracted the three-dimensional den- 
sity profiles for each component, with the parametric procedure 
discussed in Sec. 3.1 and computed the scaled central density, 
noh(z) ~ 2 , where h(zY = Q m (l + z) 3 + Qa is the ratio of the 
Hubble constant at redshift z with respect to its present value 
Hq. This parameter can be used to classify clusters into cool- 



• (noh(z ) 



> 4 x 10~ 2 ctrT 3 ) and non cool-core objects (Pratt 
|et al. 12009} . As expected, A shows a central density (noh(z)~ z = 
7 X 10" z cm" 3 ) typical of cool-core objects, while B and C show 
much lower central densities (n^hizY 2 = 2 x 10 3 crrT 3 , both). 
We extracted spectra in four (three for component C) annuli and 
fitted them with a single-temperature absorbed model, fixing as 
many components as we could because of the faintness of the 
source: nH was fixed to the galactic value ( |Dickey & Lockman 
|1990| ), redshift of the three components to 0.45 and we fixed also 
all parameters of the background components. In most cases we 
fixed also the metallicity to 0.3 Z Q , except in the centre of sub- 
cluster A where we could estimate an excess of metal abundance 
(Z = 0.6 + 0.1 Z Q ), as often found in cool cores. Due to the poor 
statistics, all temperature profiles are consistent at one <x with 
being flat and with the global values shown in Table[T] therefore 
from now on we will consider them to be isothermal. 
We combined the three-dimensional density profile and the 
global temperature to derive two other thermodynamic quanti- 
ties: pressure and entropy]^] Pressure is especially relevant to our 
analysis since it is the quantity that is measured with the SZ 
effect. We have fitted the profiles with the model described in 
Arnaud et al. (2010|l, the best fit parameters are consistent with 



the ones for relaxed cool core objects for component A, and for 
disturbed objects for component B and C. 
Entropy is a thermodynamic quantity that is connected both to 
the accretion history of the cluster and to non gravitational pro- 
cesses. If we fit the profile with a power-law plus a constant, 
the central entropy Kq is a good indicator of the cool core state 
( |Cavagnolo et al. 2009). The central entropy values are essen- 
tially driven by the central densities because we assumed a con- 
stant temperature, given the large uncertainties and poor resolu- 
tion of the temperature profiles. As expected, for subcluster A 



3 The "X-ray astronomer's entropy" is defined as K = kT/nl 13 , where 
n e is the electron density and T the X-ray temperature. This quantity is 
related to the thermodynamic entropy by a logarithm and an additive 
constant. 



5 



Planck Collaboration: The dynamical structure of PLCKG2 14.6+37.0, a Planck discovered triple system of galaxy clusters 



10.00 



10.00 t 



10. oc 



m 1-00 



S> 1.00 r* 



% 0.10 



0.10 



0.10 



0.01 



0.01 



0.01 



0." 



r [arcmin] 

A 



1 2 3 

r [arcminl 



0.' 



8 ~ 
6 - 



2 - 
L. 



r [arcmin] 



1 2 
r [arcminl 



2 - 
L. 



r [arcmin] 




1 2 3 

r [arcminl 



1000 

u 

$ 100 



10 1 

0.01 



0.10 1.00 

R/R500 



1000 



¥ 100 



10 L 

0.01 



0.10 

R/R500 



.00 



1000 



100 



1 I 

0.01 



0.10 1.00 
R/R500 



100.0 



10.0 



0.01 



.00 



R/R5 



100.0 f 



10.0 



1 .0 r 



0.1 I 

0.01 




0.10 
R/R500 



.00 



100.0 f 



10.0 



1 .0 r 



0.1 L 
0.01 



0.10 1.00 
R/R500 



Fig. 3: Radial profiles for the relevant X-ray quantities for each component. From top to bottom: surface brightness in the energy 
band 0.5 - 2.5 keV, projected temperature both as a function of the projected distance from the centre, three-dimensional entropy 
and pressure (rescaled by the value at R500) as a function of the distance from the centre in units of R500 (with the values in Table 
1). The black lines in the last two rows are the combination of the density model with temperature to estimate entropy and pressure 
and the shaded area shows the one <x uncertainty. The red and green lines are our best fit models, with the functions discussed in the 
text. 



we found ^0 = (13 + 2)keVcm 2 , a central entropy typical for 
cool core systems, while for B and C we found larger values 
(K = 142 ± lOkeVcm 2 and K = 153 ± 18keVcm 2 respec- 
tively) typical for unrelaxed objects. 



3.4. 2D structure of the components and of the supercluster 

A qualitative analysis of the X-ray image (Fig. |4]l shows that 
the two southern components are apparently connected. Indeed 
the X-ray surface brightness isophotes of component B (Fig.[T]i 
are slightly elongated in the direction of component A, as often 
observed in pairs of merging clusters (e. g. , the three systems in 
pyiaurogordato et al.|20n] and the pair A399-A401 |Sakelliou &| 
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Fig. 5: X-ray Longitudinal profile in the East- West (left panel) and North-South directions (right panel). Left panel: Negative dis- 
tances correspond to component B, positive ones to component A. The red and blue lines show the "undisturbed" models for 
components B and A (see text), while the green line in the intersection region is the sum of the two models. Right panel: Negative 
distances correspond to component A, positive ones to component C. The red and blue lines show the "undisturbed" models for 
components A and C (see text). 
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Fig. 4: XMM-Newton wavelet filtered image of 
PLCKG2 14.6+37.0. Contours overlaid correspond to the 
levels where we start to see connection between the compo- 
nents. The green and yellow regions are the ones where we 
extracted longitudinal profiles. 



Ponman 2004). We investigated visually the possible connection 
between the components by drawing constant surface brightness 
contours in the X-ray image. The appearance of the contours 
may provide information about the two-dimensional distribution 
of the intracluster gas and the possible contamination by resid- 



ual point sources. A connection between the components A and 
B is robustly detected, at a contour level above the background 
intensity (inner contour in Fig. RJ. However, with this simple 
analysis it is not possible to assess whether this connection is 
real or only a projection effect. We used the same method be- 
tween A and C where we start to see a connection at a much 
lower intensity, about 25% of the level of the background model 
in the same region (outer contour in Fig. HI. In this regime, it 
might still be possible that the connection between the two com- 
ponents is due to uncertainties in background estimation or to 
residual point sources. 

On more quantitative bases, we extracted longitudinal surface 
brightness profiles in the East- West direction across components 
A and B, and in the North-South direction, across components 
A and C (cyan and green boxes in Fig. 4}. The profile across 
components A and B (Fig. [5] left panel) shows clearly enhanced 
emission with respect to the opposite direction between the two 
clumps: we modelled the emission of each component by tak- 
ing the data in the external part of the pair and we project it 
symmetrically in the direction of the possible interaction (blue 
and red lines). In the region where the two emissions overlap, 
we summed the two models and found their sum to be consis- 
tent with the data. The two objects are very close in the plane of 
the sky and their emissions apparently overlap at less than R500 
(Fig.[T|; if they were located at the same distance from us and 
interacting we would expect to see compression and enhanced 
X-ray emission between the two objects. This is not the case 
here, so our results argue in favour of a separation along the line 
of sight of the two components, possibly still in an early phase 
of interaction. 

Concerning components A and C, their distance in the plane of 
the sky is 7.4', corresponding to ~ 2.5 Mpc at z = 0.45. The anal- 
ysis of the longitudinal surface brightness profile across them 
(Fig. [5] right panel) confirms our earlier indication: the emission 
in the intersection region is not significantly detected and is con- 
sistent with the "undisturbed" model (derived as before). 
These results suggest that the three clusters, while likely belong- 
ing to the same structure, have not started to interact yet. 
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Fig. 6: Median photometric redshifts for the three clumps (A 
black circles, B red squares, C green diamonds) of all the galax- 
ies within the cutoff radius. The error bars are the standard devi- 
ations of the photometric redshifts distribution. Given the small 
separation between components A and B, the points at radii i 2' 
may be contaminated by galaxies of the other structure. 



4.2. Optical appearance and morphology of the cluster 

We used the catalogue from SDSS to build two-dimensional 
galaxy density maps in different photometric redshift cuts (Fig. 
17), with a width &z p hot = 0.04. We assigned the galaxies to a 
fine grid of 24" per pixel, which is then degraded with a gaus- 
sian beam to an effective resolution of 3 arcminutes. We also 
computed a significance map using as reference ten random non 
overlapping control regions in a 9 deg 2 area around the system. 
Clear galaxy overdensities show up around z P hm = 0.46 at the 
location of the three X-ray clumps. However, these overdensi- 
ties do not appear isolated. At the location of cluster B, we see 
an overpopulation of galaxies towards higher redshift (5cr peak 
at Zphot ~ 0.5), consistent with the redshift z = 0.48 we found in 
X-rays, whereas the overdensity extends towards slightly lower 
redshifts (z P hot ~ 0.42) at the position of cluster A. There are 
also indications of another concentration close to component B 
at larger redshift (0.52 - 0.6). 

We investigated the maps in Fig. [TJto look for a possible popula- 
tion of inter-cluster galaxies: i. e. , objects not associated with 
one of the three clumps but rather with the whole structure, 
which would support a scenario where the three clumps are 
physically connected. We draw iso-contours levels in the sig- 
nificance map (Fig. [7): the outermost contour between 0.44 and 
0.52 connecting the three clumps indicates the presence of a 3cr 
excess in the galaxy number density above background in the 
inter-cluster region. 



4. Comparison with optical data 

Since the sky region of PLCKG214. 6+37.0 is covered by the 
Sloan Digital Sky Survejj^] we retrieved a galaxy catalogue from 
SDSS Data Release 8 (DR8). It covers a circular area of 20' ra- 
dius around the barycentre of PLCKG214. 6+37.0 and includes 
optical magnitudes and photometric redshifts (see Abazajian 
et al.|200"9 for description of measurements and calibrations of 
photometric readshifts in SDSS DR8). It contains about 2000 
objects, ~ 900 of which are in the redshift range 0.35 - 0.6. 
Unfortunately, spectroscopic redshifts are available only for the 
brightest central galaxies of components A and C, thus we relied 
on photometric redshifts alone in our analysis. Spectroscopic in- 
formation on this system will be available from our follow-up 
program and will be discussed in a forthcoming paper. 



4. 1 . Photometric redshifts of the three components 

We used the archival photometric redshifts in our catalogue to 
estimate the redshift of the three components. We extracted a 
sub-catalogue selecting only galaxies in the photometric redshift 
range 0.35 - 0.6 and, for each clump, we calculated the me- 
dian redshift of the galaxy population around the X-ray centre 
as a function of the cutoff radius. The resulting plot is shown in 
Fig. [6] Components A and C are both consistent with the spec- 
troscopic value of their central galaxies (z - 0.45), but the in- 
nermost 2' of component B indicate a slightly larger redshift 
(Zphot — 0.47), similar to the results of the X-ray analysis, al- 
though consistent at two cr with the values of the other compo- 
nents. At larger radii, the redshift estimates for the three clumps 
are all consistent with each other. It should be noted however 
that components A and B are separated only by ^ 2 arcmin, thus 
all the estimates at similar or larger radii may be contaminated 
by galaxies belonging to the other cluster. 



5. Comparison with Planck 

5. 1 . Total SZ signal 

As a simple comparison of the SZ and X-ray properties, we 
can compare the Planck Y measurement with the predicted 
values from the sum of the Y x of all three components, us- 
ing the scaling relations in Arnaud et al.| (2010). From our 
X-ray estimates (Table 1), we predict the total integrated 
value of the Comptonization parameter within a sphere of 
radius 5/?5oo for the sum of the three components to be 
Y,,5R5qq = (7.52 ± 0.9) ■ 10~ 4 arcmin 2 . This is about 50% of the 
measured signal in the same region which was found in Paper 
I and the two values are compatible at 2.3 cr. In the following, 
we will work under the assumption that the three clusters are 
all located at the same redshift z = 0.45. Considering the best 
fit redshift for component B leads to a slightly smaller SZ flux 
Y X £k5oq = (6.44 + 0.7) ■ 10~ 4 arcmin 2 . As discussed in Sec. 



2.1 we used the parameters provided in Table [T] to improve our 
estimate of the total SZ signal of this structure from Planck 
data. We built a specific template from the X-ray analysis, made 
from three universal pressure profiles cut to SR^qq (Arnaud 



http : / /www . sdss . org/ 



et al. 2010 1 corresponding to the three components. Each 
component is placed at its precise coordinates and the size is 
given by the /? 5 oo value in Table [T] We also fixed the relative 
intensity between the components to verify A/C= 0.96/1.61 
and B/C= 1.22/1.61 for the ratio of integrated fluxes. Then we 
ran the MMF3 algorithm (Melin et al. 2006) to estimate the 
amplitude of the template three times, centering the maps on 
components A, B and C. The MMF3 algorithm estimates the 
noise (instrumental and astrophysical) in a region of 10 x 10 
deg 2 around the centre (excluding the region within 5/?5oo), 
therefore changing the centering from one component to 
the other can affect the background estimation and therefore 
the flux and signal-to-noise ratio. Centering maps on com- 
ponent A we found Y 5R5m = (9.75 ± 3.19) lO^arcmin 2 , 
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Fig. 7: Galaxy density maps for cluster members in the SDSS catalogue (colours) in different photometric redshift cuts: 0.40 - 0.44 
(upper left), 0.44 - 0.48 (upper right), 0.48 - 0.52 (lower left) and 0.52 - 0.56 (lower right). The cyan contours overlaid mark the 
significance of each density peak at 3, 4, and 5 cr, respectively. The black contours show the X-ray distribution and the red circles 
and letters mark the three components. 



on B Y 5R5m = (12.26 ± 3.22) lO^arcmin 2 and on C 
Y 5Rsm = (12.97 + 3.20) lO^arcmin 2 . Our SZ flux estima- 
tions are all compatible with each other. They are also slightly 
larger than the X-ray prediction but consistent at 0.7 - 1.7 cr 
(0.9 — 1.9 cr using z = 0.516 for component E). We further 
allowed the position of the template to be a free parameter and 
found that the algorithm is able to reconstruct the position of the 
peak with a positional accuracy of one sky pixel (1.717x1.717 < 



arcrnin , in a HEALPIX projection of nside 
|2U05 i 



2048, iGorski 



et al. 



. This is consistent with the positional accuracy of 
the MMF3 algorithm, which has been tested both on simulations 
and on real data with known clusters. 

The discrepancy between the observed SZ signal and the 
prediction from the Yx measurement is decreased with respect 
to Paper I: while the X-ray prediction was only 40% of the SZ 
measurement, it is now between 60 and 77%, depending on 
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the map centering. This is partly due to the larger Fx values 
we found in this analysis with respect to Paper I, especially for 
components B and C. It is also certainly due to the improved 
accuracy of the HFI maps obtained with two full surveys of the 
sky and to the multi -component model we have used to estimate 
the SZ flux, with respect to the data from the first sky survey and 
to the single component model that was used in Paper I. Indeed, 
these results confirm our capability to extract faint SZ signals, 
when guided by X-ray priors (Planck Collaborati on|201 Te] >. 



5.2. SZ signal distribution 

It is possible to combine the X-ray images with the temperatures 
of the components to predict the distribution of the SZ signal 
(see Mroczkowski et al. 2012 for a similar approach). X-ray 



images in the soft band are proportional to the square of the 
density integrated along the line of sight and therefore their 
square root can be combined with a temperature map to derive 
a pseudo-pressure map[^] which when smoothed with the Planck 
resolution, can be qualitatively compared with the y-maps. 
We combined the background subtracted X-ray image with a 
temperature map, built assuming the mean temperature value in 
each component (Table [T| within R^m and zero outside, to pro- 
duce a pseudo-pressure map, that we smoothed with a gaussian 
filter of 10' FWHM to mimic the resolution of Planck y-maps. 
Planck cannot spatially resolve the three components of this 
object, therefore we expect the peak of the pseudo-pressure 
map to be located around the barycentre of the system just for 
resolution effects. The results are shown in Fig. [8] compared 
with the MILCA y-map. The position of the peak in the SZ map 
does not coincide with the peak of the pseudo pressure map: 
while the latter is located as expected at the barycentre between 
the three components, the y-map suggests an excess of pressure 
to the SW of component A. The offset between the two peaks is 

We have performed some tests both on the X-ray and on the SZ 
maps to investigate the origin of this offset. On the X-ray side, 
we have produced surface brightness images using a different 
background modelling. The first test concerned the background 
subtraction: we used the ESAS softwar^] to produce particle 
background and residual soft proton images and we created 
images of the "sky background" components (CXB and galactic 
foregrounds) modelling them in an external annulus ( |Leccardi| 
& Molendi 2008) and rescaling them across the field of view 
( Ettori et al.|2010| l. Point sources could also affect the position 
of the pseudo pressure peak, therefore we ran a different point 
source algorithm using the SAS task edetect_chain on MOS 
and pn images in five energy bands and we added undetected 
sources we identified with a visual inspection of the images. 
Both these tests showed a negligible impact on the position of 
the peak, which in fact is located where it is expected to be, at 
the barycentre of the three components. 

For the SZ effect, we have compared the maps reconstructed 



5 Although deriving pseudo-pressure maps as discussed in the text is 
customary in the literature, we underline here that this approach is not 
completely valid. The X-ray surface brightness, ignoring the tempera- 
ture dependance, is proportional to j n 2 dl and its square root is never 
equal to [ndl, which is the expression that should enter in the definition 
of the comptonization parameter y. However, pseudo-pressure maps can 
still be used for qualitative comparison with the y-maps. 

6 http : / /heasar c . nasa . gov/ docs/ xmm/xitimhp_xmmesas . 



html 



with different ILC-based algorithms. Besides MILCA, we tested 
GMCA ( |Bobin et al.|2008] > and NILC ( [Delabrouille et al.|2009| > 
algorithms (see Planck intermediate paper: SZ and pressure 
profile of galaxy clusters, Planck collaboration 2012, for a 
summary description and a comparison at the cluster scale of 
the three methods). The three y-maps are very consistent and 
the position of the peak does not change across the maps. 
Furthermore, we noticed that the offset was not found in the 
previous version of the maps, which was shown in Paper I. 
The presence of correlated noise in the y-map produced from 
Planck data can be a major source of error in the reconstruction 
of the position of clusters, and in particular in the case of 
low signal to noise systems such as PLCKG2 14.6+37.0. To 
quantify this error we have produced Monte Carlo simulations. 
We first estimated the noise covariance matrix on the Planck 
Comptonization parameter map of this system and we produced 
500 realisations of noise to each of which we added the expected 
thermal SZ effect from PLCKG2 14.6+37.0 to construct mock 
y-maps. Then, we estimated the position of the supercluster 
in each of these maps accounting for correlation of the noise. 
Finally, we computed the average and the standard deviation 
of the error on the reconstructed position and we obtained an 
average error of (4.5 + 2.5)'. Therefore, the offset between 
the reconstructed positions from the Planck y-map and the 
pseudo-pressure X-ray derived map of PLCKG214. 6+37.0 is 
consistent with being due to noise. The same applies for the 
separation between the peak in the Planck y- map we show here 
and the one which was shown in Paper I. 



6. Discussion and conclusion 

The first observations of the multi-wavelength follow-up cam- 
paign of PLCKG2 14.6+37.0, a triple system of galaxy clusters 
discovered by Planck, have allowed us to improve our under- 
standing of this object. With the new XMM-Newton observation 
we estimated the global properties of each component: the ICM 
temperatures range from 3.5 to 5 keV and the total masses within 
^500 are in the range 2.2 - 3 1O 14 M . We detected the iron Ka 
lines in the X-ray spectra of each component, and therefore we 
were able to confirm that components A and C are lying at the 
same redshift (z = 0.45). However, given the large angular sepa- 
ration of these two components (7.5', corresponding to 2.6 Mpc, 
in the plane of the sky), they have likely not started to interact 
yet and we did not detect significant excess X-ray emission be- 
tween these two components. For component B, we estimated 
a larger redshift from X-ray spectroscopy (z = 0.48), although 
consistent at two <x with the best fit value for component A. 
A similar indication is supported by the optical data, with the 
photometric redshifts we retrieved from SDSS DR8. However, 
given the large uncertainties of our redshift estimates (based both 
on X-rays and on photometry), a more detailed picture of the 
three-dimensional structure of PLCKG2 14.6+37.0 will be pos- 
sible only with the measurement of spectroscopic redshifts for a 
large sample of member galaxies, that is already foreseen with 
VLT in our follow-up program. 

Our redshift results are consistent with the three clusters being 
part of the same supercluster structure, that will eventually lead 
to the formation of a massive object (^ 1O 15 M ). This is sup- 
ported also by our analysis of the galaxy population with SDSS 
data: the galaxy density maps show the presence of a possible 
population of inter-cluster galaxies, significant at 3<x, connect- 
ing the whole system (Fig. [7]). However, the relaxed appearance 
of component A, its large distance (2.5 Mpc) in the plane of the 
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Fig. 8: MILCA y-map (left panel) compared with the pseudo pressure X-ray map degraded to the y-map resolution (right panel). 
X-ray contours and cyan circles indicating the three components are over-plotted to guide the eye. The green and blue circles mark 
the positions of the peaks in the y-map and in the pseudo-pressure distribution, respectively. 



sky from component C and along the line of sight from compo- 
nent B, as well as the absence of any detectable excess X-ray 
emission between the components may suggest that we are wit- 
nessing a very early phase of interaction. 

Using the X-ray results from the new XMM-Newton observation, 
we built a multicomponent model that we used to extract the to- 
tal SZ signal from Planck data. We compared the improved es- 
timate of Ysz with the prediction from X-rays and we found the 
latter to be about 68% of the measured SZ signal. The discrep- 
ancy between these two values is reduced with respect to Paper 
I and is only marginally significant at 1 2cr. 
The results from our simulations have shown that an offset 
as large as 5' can be expected in the reconstructed y-maps 
for low significance objects, due to noise fluctuations and as- 
trophysical contributions. With this study we have illustrated 
the expected difficulty of accurately reconstructing the two- 
dimensional SZ signal for objects with low signal-to-noise ra- 
tio. Indeed the instrumental noise and astrophysical contamina- 
tion compete seriously with the SZ effect at the detection limit 
threshold. Nonetheless, objects like PLCKG2 14.6+37.0 can be 
detected with a dedicated optimal filtering detection method, and 
the SZ signal can be reconstructed assuming priors (such as po- 
sition, size and relative intensity) from other wavelengths. 
Despite a deep re-observation of this system with XMM-Newton, 
the intrinsic limitations of our X-ray data and of the current 
Planck SZ maps do not allow us for the time being to assess 



the presence of possible inter-cluster emission. 
A careful analysis of the galaxy dynamics in the complex poten- 
tial of this object and of the mass distribution from weak lensing 
will both be available with our on-going optical follow up pro- 
gram. These observations, combined with the results presented 
in this paper and with new Planck data obtained in two other 
full surveys of the sky, might deliver further clues for the under- 
standing of this peculiar triple system. 
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